Quantized vortices in two dimensional solid He 
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Abstract. Diagonal and off-diagonal properties of 2D solid ''He systems doped with a 
quantized vortex have been investigated via the Shadow Path Integral Ground State method 
using the fixed-phase approach. The chosen approximate phase induces the standard Onsager- 
Feynman flow field. In this approximation the vortex acts as a static external potential and the 
resulting Hamiltonian can be treated exactly with Quantum Monte Carlo methods. The vortex 
core is found to sit in an interstitial site and a very weak relaxation of the lattice positions 
away from the vortex core position has been observed. Also other properties like Bragg peaks 
in the static structure factor or the behavior of vacancies are very little affected by the presence 
of the vortex. We have computed also the one-body density matrix in perfect and defected 
*He crystals finding that the vortex has no sensible effect on the off-diagonal long range tail 
of the density matrix. Within the assumed Onsager Feynman phase, we find that a quantized 
vortex cannot auto-sustain itself unless a condensate is already present like when dislocations 
are present. It remains to be investigated if backfiow can change this conclusion. 



Quantized vortices are one of the most genuine manifestation of the presence of superfluidity 
in many body quantum systems, but from the microscopic point of view no complete 
understanding of them has been reached yet. Recently quantized vortices have been related 
to the supersolidity issue [Jj; in fact, arguments in favor of a vortex phase in low temperature 
solid ^He, preceding the supersolid transition, are appeared in literature [2l[3lll]. With respect 
to this possible connection one of the fundamental questions to be answered is: what does a 
quantum vortex look like in solid Helium from a microscopic point of view? 

Dealing with vortices is a really hard task for microscopic methods, and it calls for some 
approximations or assumptions [5-10]. In fact, the wave function has to be an eigenstate of the 
angular momentum, so it needs a phase. Following the well established routine for the ground 
state, once chosen a variational ansatz for the wave function, one could be tempted to correct it 
by means of exact zero temperature Quantum Monte Carlo (QMC) techniques. Unfortunately 
this is actually not possible because of the sign problem that affects QMC methods. The most 
followed recipe is to improve the variational description via QMC, but releasing the exactness 
of the methods in favor of approximations that allow to avoid the sign problem, like for example 
fixed phase [6] or fixed nodes [10] . 

We study here the properties of a single vortex in solid ^He via the Shadow-PIGS (SPIGS) 
method with fixed phase approximation. The many-body wave function can be written as 
^(R) = e^^(^)^'o(i2), where n{R) is a many-body phase, ^o{R) is the modulus of the wave 
function and R = (ri, r2, . . . , r/v) are the coordinates of the N particles. ^{R) describes a 
quantum state of the system if it is a solution of the time independent Schrodinger equation: 
from H'^{R) = E^{R) it is possible to obtain two coupled differential equations for and 
for ^'o(i?). The fixed phase approximation consists in assuming the functional form of Q{R) as 



given and to solve the equation 



=1 1=1 



(1) 



for "i^Q^R). Solving ([T]) is equivalent to solve the original time independent Schrodinger equation 



for the A^-particle with an extra potential term = 2m ^i=ii^i^) ■ 

The simplest choice for the phase is the well known Onsager-Feynman (OF) phase [11]: 
= lJ2iLi(^i (where 9i is the angular polar coordinate of the i-th particle). ^{R) is an 
eigenstate of the z component of the angular momentum operator with eigenvalue hNl, 
being / = 1,2, .. . the quantum of circulation. This choice for 0,{R) gives rise to the standard 
OF flow field: in fact the extra-potential in ([1]) reads 
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where is the radial polar coordinate of the i-th particle. 

In order to sustain a quantized vortex, the system should display a macroscopic phase 
coherence, and at T = K this means that solid ^He should house a Bose-Einstein condensate 
(BEC). It is known from QMC results that no BEG is present in the perfect crystal [12-15], but 
if the vortex turns out to be able to induce a BEC it could be a self-sustaining excitation. On the 
other hand, it is largely accepted that defects are able to induce BEC [16j, and then a defected 
crystal can safely sustain a quantized vortex. Here we report on the study of a two dimensional 
(2D) ^He crystal with and without dislocations. In fact, dislocations can be included in the 2D 
crystal without imposing boundary constraints [16]. Moreover the 2D system allows to reach 
large distances keeping the number of particle in the simulation at a tractable level, and this is 
a desirable feature when interested in off-diagonal properties of the system. 

We face the task of solving ([T]) with the extra-potential given by ^ when I = 1 with the 
SPIGS method \17\ [TH] , which allows to obtain the lowest eigenstate of a given Hamiltonian 
by projecting in imaginary time a SWF [TO] taken as trial wave function. The SPIGS method 
is unbiased by the choice of the trial wave function and the only inputs are the interparticle 
potential and the approximation for the imaginary time propagator [20j. As He-He interatomic 
potential we have considered the HFDHE2 Aziz potential [21] and we have employed the pair- 
Suzuki approximation [20] for the imaginary time propagator with time step St = 1/120 K^^. 
One difficulty with ([2]) is that the potential is long range so that either one puts the system in 
a bucket [6l [7] or one should consider a vortex lattice ^ . Such complications can be avoided by 
multiplying l/r^ in ([2|) by a smoothing function 
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(L being the side of the simulation box) so that standard periodic boundary conditions can 
be applied. With this choice, the extra-potential is equivalent to the OF one only for r < A, 
the provided ^'o is no more an exact eigenstate of Lz but it is close to it in the interesting 
region of the vortex core if A is large enough. Here we have used A = 30A. We have performed 
simulations at p = 0.0765A~^ in a nearly squared box designed to house a perfect triangular 
crystal with M = 572 lattice sites, and a crystal with 10 vacancies (A^ = 562); such vacancies 
in the initial configuration transform themselves in dislocations |16j . 

In Fig. [T]we report our results for the integrated vortex energies et,(r) = (Ey{r) — E(r))/N 
{E^{r) and E{r) are, respectively, the energy of the particles that lie inside the disk of radius r 
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Figure 1. Integrated vortex energy 
e„ with error bars as a function of the 
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Figure 2. Comparison between the radial 
defect distribution pd{r) as a function of 
the distance from the origin (vortex core 
when vortex is present) for a 2D sohd '^He 
si p = 0.0765A-2 with (fihed symbols) 
and without (open symbols) a vortex 
in the perfect (circles) and in defected 
(squares) crystal. 



in the system with and without vortex) as a function of the distance from the core in the perfect 
crystal. The center of mass of the system is not fixed, and we find that, independently on the 
starting configuration of the crystal, the vortex core sits in an interstitial site. We also find a 
very small relaxation of the surrounding lattice around the vortex core. 

In order to study the effects of the vortex on the crystal properties we have monitored the 
static structure factor, the pair distribution function and the radial defect distribution pdir) 
(i.e. the distribution of the particle whose coordination is different from 6 as a function of the 
distance from the origin where the vortex core is located). In Fig. [2] we plot our results for pd(r) 
both for the perfect and for the defected crystal. We find that in the defected case the radial 
defect distribution is about an order of magnitude larger than in the perfect one; however, the 
results of the system with and without vortex are very close each other so that we conclude that 
the OF vortex does not affect in a sensible way the disorder which behaves as in the system 
without vortex. We come at a similar conclusion for the crystalline structure, since both the 
static structure factor and the pair distribution function show no appreciable differences for the 
system with and without the vortex. 

We have computed also the one-body density matrix pi for the perfect and for the defected 
crystal with and without the vortex in order to investigate the vortex effects on the off-diagonal 
properties. Our results are reported in Fig. [3j pi for the crystal with vortex are indistinguishable 
within the error bars from the ones obtained without vortex. We conclude that the OF vortex 
is not able to induce a Bose-Einstein condensation (BEC) in the perfect crystal, or to increase 
the already present condensate fraction in the defected one. 
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Figure 3. One body density matrix pi 
in 2D solid ^He without (open symbols) 
and with (filled symbols) a vortex for the 
perfect crystal (circles) and for the crystal 
with dislocations (squares). 



Since no BEC is present in the perfect crystal and the OF vortex is not able to induce it by 
itself, we can conclude that perfect 2D solid ^He can not sustain vortices of the OF type. Thus 
the OF wave function is a possible representation of a vortex in solid ^Helium only when BEC 
is already present, like in a defected crystal. Preliminary results in three dimensional solid ^He 
seem to confirm such conclusions for the OF vortex. 

In the liquid phase, the OF phase has been improved with the inclusion of back-flow 
(BE) correlations [6l [8]. The effect of backflow increases [8j at higher densities and it might 
well become dominant in the solid phase. Computations along this line are in progress. Since 
the BE terms acts mainly near the vortex core, we might expect that BE will not modify much 
the off-diagonal properties in 2D, where the vortex core is a point defect, while it could become 
relevant in 3D where the vortex core is an extended defect. 
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